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We consider a system of interacting particles subjected to Langevin inertial dynamics and derive 
the governing time-dependent equation for the one-body density. We show that, after suitable trun- 
cations of the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy, and a multiple time scale analysis, 
we obtain a self-consistent equation involving only the one-body density. This study extends to ar- 
bitrary dimensions previous work on a one-dimensional fluid and highlights the subtelties of kinetic 
theory in the derivation of dynamical density functional theory. 
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I. INTRODUCTION 

Equilibrium density functional theory (DFT) provides a very powerful tool to investigate the structure of inhomo- 
geneous fluidsis^!^. For systems at thermodynamic equilibrium, all the statistical properties can be expressed in terms 
of the average density field, p(x). In turn, p(x) can be obtained by minimizing the grand potential functional, f7[p(x)], 
which for a given interaction potential and external one-body potential, is unique. Although ri[p(x)] is known exactly 
only in few particular cases, fairly good approximations can be devised, so that the method is versatile and generally 
applicable with success. 

In out of equilibrium situations, there is no rigorous analogue of the free energy, so that the description must be 
necessarily more complex. In many instances, the density field is not sufficient to represent the state of a system, 
and one has to identify the minimal set of relevant fields which provide a satisfactory description of its evolution. 
In the past, phenonienological equations have been devised to deal with a large class of non equilibrium processes, 
such as Navier-Stokes, reaction-diffusion, Ginzburg-Landau equations, eto^i^i^. Standard hydrodynamics, for instance, 
identifies d + 2 relevant fields, number, momentum and energy density, where d is the spatial dimensionality of the 
embedding space, each of these fields satisfying a conservation law, and thus relaxing slowly toward equilibrium. On 
the other hand, kinetic theory provides a truly microscopic derivation of the Navier-Stokes equations and requires 
some assumptions about the evolution of the phase space distribution^. 

When the dynamics of the particles is dissipative and overdamped, as in the case of colloids strongly coupled to 
the solvent, one may assume that the evolution of the density field is sufhcient to capture the relevant physics of the 
system. This occurs because the momentum and energy fluxes are not conserved any more and thus become rapidly 
enslaved by the density field. If the friction coefficient, 7, is large, the particles velocities rapidly relax in a time 
of order 7""'^, so that the evolution on longer timescales is only determined by the dynamics of the configurational 
contribution to the distribution. This is the point of view of the Smoluchowski description^iiSiiS, which takes into 
account only the configurational degrees of freedom. Based on these considerations, in recent years a dynamical 
density functional theory (DDFT) has been put forward by several authorsi^J^. In DDFT the evolution of the noise- 
averaged, time-dependent, one-particle density is completely deterministic and is driven solely by the gradient of the 
equilibrium free energy functional. 

In more general terms, inertial effects can play a role in a number of different physical systems, as in the case 
of colloidal particles weakly coupled to the solvent. Inertial effect must be accounted for by means of the Kramers 
description in full phase space^^i^i. In the present paper, we derive a deterministic evolution equation for the density 
field starting from the Kramers equation. A second goal of the present derivation is to show that the kinetic approach 
leads quite naturally to an evolution equation which contains DDFT as a limiting case. 

There are a few reasons for working in the framework of the Kramers equation. The first is that Kramers is a more 
general description since it contains as limiting cases the Liouville equation and the fully overdamped dynamics. In 
addition, the Kramers equation allows to deal quite naturally with systems under temperature gradients^. Finally, 
from a fundamental point of view, it seems to us rewarding to reconcile the density functional approach, whose major 
success has been in explaining equilibrium phenomena, with the kinetic approach which is usually applied to non 
equilibrium situations. Whereas the current treatments show that it is quite natural to obtain the DDFT equations 
for overdamped dynamics, it is not straightforward to derive the DDFT from the Liouville equation, such as in the 
case of a dense liquid, where the hydrodynamic currents carried by particles are not strongly damped by the viscosity 
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of the surrounding fluid. A previous work by Marconi and Tarazona32i has shown that in the case of one dimensional 
hard-rods the full inertial dynamics gives rise to DDFT in a perturbative analysis truncated at lower order. However, 
the extension to d-dimensions of this procedure is not trivial, because the number of independent tensorial components 
of the moments of the distribution increases with the order, n, as d". 

We briefly recall the several approaches appeared in recent literature to build the equation of evolution of the one- 
body density profile in a classical fluid. One class of approaches starts from the fundamental Hamiltonian description 
of the fluid and aims to derive an equivalent evolution equation for the density2i. By using the Kawasaki-Gunton 
projection operator formalism^^,^ Yoshimori^? has recently derived an equation of motion for the density which differs 
from standard DDFT and is stochastic on account of the coupling with fast degrees of freedom not explicitly retained. 
Another class of equations for the density is obtained starting from the Langevin microscopic dynamics of the particles. 
Within this class, some authors have kept the stochastic character of the evolution of the instantaneous densit y^^ 
whereas others have focused attention on an ensemble averaged density and obtained a purely deterministic equation. 
Both methods do not require spatio-temporal coarse graining^^. 

The different approaches underline the need for a theoretical framework which is capable of deriving equations 
similar to the DDFT in systems described by second order thermostatted inertial dynamics, instead of the Langevin 
overdamped dynamics. To this purpose one must consider a phase space rather than a configurational space de- 
scription. The natural framework turns out to be the kinetic theory which has predated the modern liquid state 
approaches. Unfortunately, a Boltzmann equation suitably modified to include the coupling to the heat bath only 
applies to a dilute gas regime. However, several authors^^i^i^ have extended the kinetic approach to the dense 
regime, where the harsh repulsion between particles is taken into account via a revised version of the Enskog theory 
(RET). Ernst and van Beijeren have obtained the RET equation for the distribution function from the full phase 
space representation of the N-particle dynamics. In addition, Ernst and van Noije have considered the evolution of 
thermostatted systems-^. We assume this theory as the starting point to derive the time evolution equation for the 
one body density for hard-sphere systems. 

The present paper is organized as follows: in section |lT] we describe the microscopic model and write the evolu- 
tion equation for the N-particle distribution function. Subsequently, we contract the description to the one-particle 
distribution function, using the BBGKY hierarchy^, and separate the velocity from the spatial dependence of the 
distribution obtaining an open hierarchy for the moments. In section IIIII we close the hierarchy of the distribution 
functions in the case of two explicit models for the intermolecular forces. In the case of soft potentials we apply a 
Random Phase Approximation, whereas in the case of hard-core potentials we consider the Revised Enskog Approxi- 
mation^!^. In section ITVl we apply the multi-scale method to close the open hierarchy of the moments and to obtain 
an equation for the density of the system. 



II. MODEL 



Let us consider an assembly of N heavy particles suspended in a solution of light particles and moving in a d- 
dimensional region with positions r„ and velocities u„ with n = I, N. Due to their small mass, the solvent particles 
perform rapid motions so that their influence on the heavy particles can be described according to Langevin's idea 
by means of an effective stochastic force. As a result of such elimination of microscopic degrees of freedom the heavy 
particles experience a viscous drag force proportional to their velocity plus a random stochastic acceleration. In 
addition the N heavy particles of identical mass m may move under the action of an external force fe(x) and interact 
through a pair potential U{\r — r'|), which for simplicity we take as a continuous function'^^-^^. The equations of 
motions are: 



^dt 



fe(r„) - ^ Vr„C/(|r„ - r,„|) - TO7U„ 



dt + dWn{t) 



(1) 
(2) 



where (iW„(i) = ^„(t)di is the increment of the Wiener process and ^„(i) is a Gaussian white noise with properties 



(C(i)> = 



(3) 



where T is the "heat-bath temperature" and (•) indicates the average over a statistical ensemble of noise realizations^. 
We include the friction exerted by the solvent via the coefficient 7. 

We now derive an equation for the singly conditioned probability distribution p({x, v}, t|{xo, vq}, 0) defined as the 



probability that the random variables {r{t),u{t)} will lie between {x, v} and {x+dx, v + dv} given that {r(0), u(0)} = 
{xo,vo}^. 

Let us consider the infinitesimal time evolution of an arbitrary function G{{r{t),u{t)}) of the trajectory of the 
system^ 



dG({r(i),u(i)}) = 5]{[u„.Vr 



y Vr„C/(|r„-r^|)-7U„)-Vu„ 
Vu„ • V„ J dt + -dW„(t) • Vu„ }G{{r{t),u{t)}) 



The average over the realizations of such a function reads 

d 



dt 



{G{{r{t),u{t)})) 



n 



n ' ^ r. 



'fe(r„) 1 



TO m 



Vr„C/(|r„ -r^l) -7U„yV 



G{{r{t),n{t)}) 



(4) 



(5) 



Since the random variables {r, u} have the conditional probability density p({x, v}, i|{xo, vq}, 0), we can equally 
write 



(G({r(i),u(t)})> - / d^^xd^%G({x,v}M{x,v},t|{xo,vo},0) 



(6) 



and 



|(G({r(t),u(t)})) = I d^xd^vG({x,v})|p({x,v},i|{xo,vo},0) 
Similarly, the r.h.s. of equation ^ can be expressed as 

J d^xd^v^l [v„ . V.„ + - E Vx„t/(|x„ - x„,|) - 7v„)-V. 



(7) 



TO 



Vv„ • Vv„] }g({x, v})p({x, v}, t|{xo, vol, 0) 



(8) 



We next integrate by parts and discard surface terms in equation fS]) and use the arbitrariness of the function G. The 
result is 

d 



- + E[4"^-4'p])M{x,v},t|{xo,vo},0) 



Vx„?7(|x„-x™|)- Vv„p({x,v},t|{xo,vo},0) 



where we have introduced the notation 



r(") _ V V 



fe(x„) 



r(n) 



(9) 

(10) 
(11) 



To make contact with the standard statistical mechanical approach, we introduce the F phase-space distribution, 
/(^)({x, v}, t) as the probability density of finding the many body system between {x, v} and {x + dx, v + dv} at 
time t. Knowing the value of f''^^ at time < = we obtain its value at time t via the relation 



/W({x,v},t)- / dxodvop({x,v},<|{xo,vo},0)/W({xo,vo},0) 



(12) 
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One concludes that the evolution equation for /*^^)({x, v}, i) has the same form as equation ^ relative to 
p({x,v},i|{xo,vo},0), i.e.: 

(I + E[4"^-4t])/^"^({x,v},i) = 

71 

-EE Vx„{/(|x„-x„|).Vv„/W({x,v},t) (13) 



TO 



If one sets 7 = 0, the term L^pp drops and the dynamics turns out to be conservative and the evolution of /^^^ is 
described by Liouville equation. 

The description contained in /^^^ is for any practical purpose unnecessarily detailed, so that we are lead to consider 
the reduced s-particle distribution functions obtained by integrating over 2d{N — s) degrees of freedom: 



/('*)(xi, .. , Xs,vi, ..,Vs,i) = 



{N ^s)\ 



n=s+l 



^=s+l 

/(^^(Xi, ..,Xs,Xs+i..,XAr,Vi, ..,Vs,Vs+i..,VAr,i). (14) 

Thus integrating equation l|13p over (N-1) particle's coordinates and velocities we obtain the exact evolution equation 
for the one particle distribution 

-/(I) (x, V, t) + L^f^^^ (x, V, t) = Lppf'^ (x, V, t) + fc(x, V, t) (15) 

where, the left hand side of equation lfT5|l contains the free streaming of the particles, while the right hand side 
describes the interactions with the heat bath and those among the particles. The interaction term fc(x, v, t), in the 
case of continuous potentials, reads 

A:(x,v,i) = ^Vv^ d^ j dv7(2)(x,v,x',v',t)VxC/(|x-x'|) (16) 

For particles with hard core interactions the term VxJ7(|x — x'|) is undefined and we have to modify the treatment 
as briefly explained in section Hill We remark that /^^^ is the p.d.f. in the 2c? dimensional /i-space and equation (jisp . 
supplemented by equation ([T6|) . represents the first equation of the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) 
hierarch}*^, which connects the evolution of the n-particle distribution function, to the distribution function for (rt + 1) 
particles. Since /^^^(x,v,f) depends on the two-particle distribution /(^^(x, v, x', v', i), some approximate closure is 
required to obtain a workable scheme. 

Before closing this section three remarks are in order: i) the stationary solution /^^^ of equation (fT5)) reduces to 
the equilibrium Gibbs (N,V,T) phase-space distribution, ii) if we take the limit 7 = 0, i.e. decouple the system from 
the heat bath, the dynamics becomes micro-canonical, i.e. at constant energy, and the evolution of Z*-^^ is completely 
reversible, iii) The number density p(x, t) can be obtained from the ensemble average of the microscopic operator 
X]„<5(x — Yn{t)), so that p{'^,t) = / dv/(^)(x, V, i), the average particle current is j(x,f) = / dv/(^^(x, v, t)v and so 
on. 

In the following we shall employ the non dimensional set of variables which arc obtained by measuring the velocities 
in units of the thermal velocity vt = y/kpT /ra and lengths in unit of ct, i.e. V = and X = ^. The remaining 
variables can be non-dimensionalized according to the transformations t = t— , F = 7—, F(X) = — 2^fe(x). 
Finally, the distribution function and the collision term are rescaled according to the transformations: P(X, V, r) = 
ct''w^/(-^)(x, V, t) and A'(X, V,t) = a'^^^Vp^^k{x,v,t), where d is the dimensionality of the embedding space. 

By rewriting equation (|15p as 



1 aP(X,V,r) 
F dr 



CppPiX, V, r) - - V • VxP(X, V, r) 
iF(X,T) • VvP(X, V,t) + iif(X, V,t) 



(17) 
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one sees that, when the non dimensional damping constant F is large, the Fokker-Planck term 

£j.pP(X, V, r) = Vv ■ [Vv + V] P(X, V, r) (18) 

is dominant. This operator has non positive eigenvalues = 0, —1, —2, .. and corresponding eigcnfunctions H^^V) 
defined as 

HlHV)^hlHVMV) (19) 

where 

= (^^"'"^ 

and the subscript a represents an array of i/ integer elements each ranging from 1 to d^^. The properties and the 
explicit representation of some of these polynomials are given in appendix A. 

To separate the spatial dependence from the velocity dependence we expand the distribution function as : 

P(x,v,t) = $(")(x,r)ffW(v) + $f^(x,r)iff)(v) 

+ l<i>g)(x,r)ij(;)(v) + ici>g) (x,r)ijg(v) + ••. (21) 
where the Einstein convention on repeated indices is used. Similarly, we expand the interaction term and obtain 

i^(x,v,r) = )(x,r)i/ii)(v) + Ic^;' (x, r)i/f (v) + l^cgl(:^,r)Hl^l(y) + ... (22) 
The first term C(°) = since the number of particles is conserved, while the remaining are obtained from the formula: 

Cf(x,r)^ y'dV/i(^)(v)if(x,v,T). (23) 

Using the orthogonality of the tensorial Hermite polynomials one derives a set of differential equations for their 
coefficients $^''(X, r). This is achieved by multiplying both sides of equation (fl7|) by and integrating with respect 
to V. One finds the following set of equations, which reduce to the standard hydrodynamic ones for F = 0: 

(X,r) = -dxM'\^.r) (24) 



i5x,(<i>g'(X,r) + ci>(f (X,r)) +cP^(X,r) (25) 



^^-<i>p)(X,T) = -F<i>f)(X,r)- A$("HX,r) 



^$(f(X,r) ^ -2F<J>(f(X,r)-[AM')(X,r)+D,$«(X,r) 



ax.$gl(X,r)+c{;)(X,r) (26) 



having introduced the short hand notation: 



III. EXPLICIT MODELS 



In order to proceed further we need to provide an explicit representation of the collision operator A'(X, V,t). We 
shall do that in two cases which may serve to illustrate different aspects of the dynamics. 
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A. Mean-field model 



Given a continuous and differcntiable pair potential J7(|X — X'|), we can write the collision kernel as 

K{X,V,t) =Vv ■ J dV' j dX'F(2)(x,X',V,V',r)Vx;7(|X-X'|) (27) 

which contains the two particle distribution function P*^^' (X, X', V, V', t). We assume, as Enskog did, 
that the doublet distribution function P^^) (X, X', V, V, r) factors into the product: P^^) (X, X', V, V, r) = 
P(X, V, t)P(X', V, t)(72(X, X'Ip) in which the pair correlation 32 is not a function of the velocities, but only of 
the particle positions. It is assumed that 52 is a nonlocal equilibrium function of the profile p(X, t) and depends on 
time only through the density profile, and to have the same form as in a nonuniform equilibrium state whose density 
is p(X, i). Given the above, 52 can be determined exactly from the second functional derivative of the free-energy 
functional for a inhomogeneous systen*^. 
With this approximation we find: 

If (X, V, t) = -F™°' (X, t) • VvP(X, V, r) (28) 
where, using p(K,T)a'^ — J (iV'P(X', V', r), we have introduced the molecular field 

F™°'(X, r) = -J dX'p(X', r) 52(X, X'jp) Vxf/(|X - X'|) (29) 

In order to obtain the coefficients Ca^ one must take the scalar product of the operator given by equation (|28|) with 
/i2^(V). One finds: 

C(°)(X,t) = (30) 
Cf)(X,T) = -$W(X,r)P™°'(X,r) (31) 
(X, r) = - ' (X, r)P™°' (X, r) + <i>f^ (X, r)Fr' (X, r)] (32) 

Let us show that at equilibrium the molecular field can be expressed as a gradient. Following Marconi and Tara- 
zonai^, we consider the relation, valid at thermodynamic equilibrium: 



/ 



dX'p{X',T)g2{X,X'\p)dx,U{X-X')=dx ' 



6piX,T) 

(33) 



where J^'^^'^ is the free energy excess over the ideal gas, so that the molecular field is of the form sought. Consequently, 
we express the first non vanishing component of the interaction as 

CP(X,r)=.VX,r)ax.^^^»^ (34) 

dp{X,T) 



B. Hard sphere model 

The hard sphere (HS) model can be viewed as a collection of hard, perfectly elastic objects. There is no inter- 
molecular force acting on them except at the instant when two spheres are at distance cr, the HS diameter. Thus HS, 
due to the coupling with the stochastic heat bath, behave like non-interacting Brownian particles between collisions. 
When two HS collide their velocities after the impact, denoted with a prime, are related to those before the impact 
(unprimed) velocities by 

V'l = Vi - (Vi2 • fi)fi 

Vj, = V2 + (Vi2 • fi)fi (35) 

where V12 = vi — V2, fi is the unit vector directed from particle 1 to particle 2. 
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Using the fact that, although the derivative of the intermolecular potential is ill-defined, the trajectories in phase 
space are well defined, Ernst and van Beijeren were able to derive the evolution equation of the N-particle distribution 
function They proved that the evolution of the T phase space distribution for a randomly driven HS fluid is 

governed by the following pseudo-Liouville equation 

n n n<m 

(36) 

The second term on the right hand side of equation (j36[) contains the binary collision operator Timn) for HS which 
accounts for the impulsive hard core interactions and reads: 

T(mn) = a'^-i j di3e(v,„„6)(v,„„ • i3) [^(x„.„ - fi)6;:„ - 5(x„,„ + fi)] (37) 

where is the Heaviside function which selects only particles approaching one another, v^n = (vm — Vn), Xmn = 
) and 6*^„ is the scattering operator defined for arbitrary function 5(v„i,v„) by 

C„5(v™,v„) =5(v;„,v;) (38) 



which for elastic HS acts on the velocities v„, and replaces them by the restituting velocities as defined by equation (j35p . 

Again starting from equation ([55]) one obtains the first level of the BBGKY hierarchy for the reduced one-particle 
distribution function /'^^(x, v,i) which is similar to equation (|15p . In the HS case the derivation is quite lengthy so 
that we merely quote the final result of Ernst and van Beijeren for the interaction term 

fc(xi,vi,i) = a'^-^ j dw2 y"d6e(fi-vi2)(B-vi2) X 

{/(2)(xi,v'i,Xi -fi,V^,t) - /(2)(xi,Vi,Xi +fi,V2,t)} (39) 



where the primes on the velocities denote scattered values determined from equation p5|) . Finally, we consider the 
so-called Revised Enskog Theory (RET) which consists in replacing Z'^', using the molecular chaos hypothesis^ 

/'^^(Xi, Vi,X2, V2,i)(5(|xi2| - 0-) = /^^^(Xi, Vi,i)/(^)(X2, V2,t)52(xi,X2|p)(5(|xi2| " cr) 

where (;2(xi, X2I/3) is the pair correlation function. It is worth noticing the difference with the mean- field approxima- 
tion. Here, it is assumed that atoms are uncorrelated immmediately prior to collision but are correlated after they 
collide^, because collision itself generates correlations. 

We turn hereafter to the nondimensional notation and use capital letters for the reduced variables. With the 
substitution (|40|) the collision integral iir(X, V, r) can be represented by the following convolution product: 

i^(Xi,Vi,T)= j dY2 y"dfie(6- Vi2)(6- V12) X 

{.92(Xi, Xi - fi|p)P(Xi, V'l, t)P(Xi - 6, V^, t) - 

.g2(Xi,Xi+fi|p)P(Xi,Vi,T)P(Xi +fi,V2,r)} (40) 

We now turn to the evaluation of the coefhcients of the collision integral via the Hermite expansion (|23p . Since this 
expansion involves an infinite number of coefficients, we need to resort to a suitable truncation scheme. A physically 
motivated prescription consists in choosing a restricted basis in the complete Hilbert-space of the functions h''a ■ To 
this purpose we choose the following 2 + d functions 

{/i(°)(v),;iW(v),/i(2)(v)} (41) 



with h^"^^ {y) = (V), corresponding to the 2 + d hydrodynamic modes (density, current and energy) which in 



4?( 

the absence of thermostat are conserved. 

In order to project i^(X,V, r) on this basis it is convenient to replace (Vi,V2,f5) — ^ (V'j^, V2, — fS) so that the 



8 



projection can be written as 



Cf(X,r) ^ J dViJ dV2 J d6e(fi • Vi2)(fi • V12) x 

[hl^ (V'l) - (Vi)]52 (X, X + 6|p)P(X, Vi , t)F(X + 6, V2, r) 



(42) 



Equation (|42|) still represents a cumbersome (Sd— l)-fold integral. However, by a change of variables and using simple 
properties of Gaussian integrals it is possible to perform 2d integrations. In order to see that, we introduce the relative 
velocity and the center of mass velocity: 



Z = Vi - V2 
Y ^ V1+V2 



(43) 



and notice that the Jacobian of transformation from (Vi, V2) to (Y, Z) is unity, so that the integral becomes, using 
the Hermite representation of P(X, V, r). 



c(')(x,- 



d%2(X,X + fi|p) / d'^Y-^e-^' I d'^Z 



1 



e 4 



e(fi- z)(6- z) 



/^£Hy + |-(6.Z)6)-/»W(Y + |) 



$ W (X + fi, r) + y hf^ (Y - !)$(>' (X + 6, r) + h^^^ (Y - (X + fi, r) 



$W(X,r) +y /.«(Y + |)$«(X,r) + (Y + |)|.(^) (X, r) 



with = 'i'L^''- The integration over Z and Y can be performed straightforwardly with the results: 

C(")(X,r) =0 



(44) 



(45) 



~ J d'^<T<7,52(X, X + fi|p){$(")(X, r)$(0)(X + fi, r) 



cW(X,r 
2 



i=i 



(X, r)$(^) (X + 6, r) + l-^^) (X, t)$(") (X + fi, t 



E 



^<5,,$f ) (X, r)$« (X + 6, r) - ^a,a,<i>« (X + fi, r)<i>« (X, r) 



1 

.7=1 



(46) 
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C(2)(X,t) = d'^ag2(X,X + 6|p) x 

d 

{ [*?^ ^)^^°^ + 6, r) + *f ^ (X + fi, t)$(°) (X, , 

^ (X, r)$(2) (X + 6, r) - l>(2) (X, t)$(°) (X + 6, t)] + 

d 

Y,&, [$f)(X,T)l>(2)(X + 6,T) + $f^(X + fi,r)|.(2)(X,r)]} (47) 
i=i 

Having now obtained the expressions of the coUision matrix elements in the restricted hydrodynamic space, we could 
in principle insert these into the moment hierarchy equations (1241126^ . However, the hierarchy is not closed and we 
need some prescription in order to truncate the expansion. In the following we shall consider a perturbative method 
to achieve this goal. 

IV. MULTISCALE METHOD 

The multiple time-scale approacb^^^ represents a powerful perturbative method to determine the temporal evolu- 
tion of the distribution function P(X, V,t) in the high friction limit, F^^ << 1. Physically, such a limit means that 
many solute-solvent collisions take place while a solute particle travels a distance of the order of its radius. The method 
is perturbative and the relevant expansion parameter is the inverse dimcnsionless friction, T^^ . By means of such 
a procedure Titulaer— , starting from the Kramers equation, has derived the Smoluchowski equation for a colloidal 
particle. The latter equation could have also be obtained, without resorting to the multiple time scale method, by 
simply taking the high friction limit in the microscopic Langevin equation, i.e. dropping the inertial term. However, 
the perturbative method allows to compute the corrections to the Smoluchowski equation when F is finite^. 

To begin with, the moments , their time derivatives and the collision integral are expanded in powers of 
F-i « 1: 



d _ d Id Id 



(48) 



= + (49) 

Let us remark that, in order to construct the solution, we have replaced the single physical time scale, r, by a 
set of auxiliary time scales (tq, ti, .., r„) which are related to the original variable by the relations t„ = F~"r and 
are treated as independent variables. Also, all time-dependent functions of r are replaced by auxiliary functions of 
(to,ti,..). Once the equations corresponding to the various orders have been determined, we return to the original 
time variable and functions. 

Following the method of reference^, the amplitudes (j)n^ with n > are set equal to zero. Such a choice, although 
not unique, is sufhcient to eliminate secular terms, i.e. terms containing a dependence on the slow time tq. 

Expansions (|48ll50p must, now, be inserted into equation ([T7| in order to identify on both sides of the resulting 
equation terms belonging to the same order in the expansion parameter F~^ and in the Hermite polynomials. We 
do not report the lengthy but straightforward derivation but merely write down the set of equations which allows to 
relate all the amplitudes (l>n}a in terms of the evolution of the amplitude ^q"^ . 

The order F^ yields: (/)^°^ ^ and cj)'^^ = > 0. 

The first power F~^ instead gives three equations, which allow to determine the amplitudes of order 1: 



(51) 



10 



and 



^ - D.4>':' - eg (52) 



= (53) 



derivative of the amplitude 0^^^ 



The following order F ^ gives and equation for the partial ti derivative of the amplitude 0q°^ in terms of a spatial 



Thus, using equation (|52p . we obtain 



At the same order we also find the relation 



5ro 



1,. - -^2. - <^il = (56) 



which using the fact that 0^^' does not depend on tq, being a functional of (/)g°'', in view of equation (|52p leads to the 
equation: 

4^ - (57) 

Finally, we find: 

= V/ > 1 (58) 

To order F"'^ we need to consider only the relation: 

^/^^ ^ -d.A'} - -dxA'J (59) 

which allows to express the partial derivative of ^o"-* with respect to T2 in terms of the collision integral. 
We restore, now, the original time derivative: 

9r'^° "^F9ri F2aT2^'^o ^^"^ 



and using equations ()55|) and (I59p we arrive at 



^0r(X,r)4vx 



( Vx - F(X))0W (X, r) - cW (X, r) - lc« (X, r) 



(61) 



As one can see, at this order in the perturbative expansion only the component = 1 of the collision operator appears 
in equation (j6ip . We need to give an explicit representation of in order to have a closed evolution equation. 
In the Vlasov model the procedure is straightforward: in fact, given eq. (|3ip. the amplitudes c^^] only depend on 

the (/)[)*''' , i.e. the density profile. Since such a quantity is of order F*^ in our small parameter expansion, it follows that 
there is no correction of order F~^. Only at order F"^, which is beyond the scope of the present treatment, one finds 
a correction to the evolution equation in the Vlasov model. 
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In the hard-sphere model, instead, we have the two following terms: 

4^](X,t) = - y dfia,g2(X,X + i3|0[,°))4"Hx,r)0(°)(X + J3,r) (62) 

c^^}{X,t) = --^ J dfia,52(X,X + 6|4°') X 
d 

a, [^15 (X, r)^'') (X + fi, r) - cjQ (X + fi, r)^") (X, r)] 

(63) 

with 

0W(X,t) = -(ax. - ^^.(X))4'^)(X,r) +c(^,)(X,t). (64) 

We show next that if we keep only the term Cq\- and drop the term c^^^l of order in equation (j6ip we recover 
the time dependent DDF equation. In fact, using the following expression for the HS impulsive force: 

ax,f/Hs(|X - X'l) = -a,S{\X - X'l - 1), 

where ?7//s(|X — X'|) is the HS pair potential, we can rewrite eq. ((62|) as 

4^](X,t) = J dX'g2(X,X'|4°))4")(X,T)4°)(X',T)ax,C///s(|X- X'l) (65) 

or, using relation (|33|) . in the equivalent form: 

c(:)(X,r).0r(X,r)ax.^gM (66) 

(50^, '(X,r) 

where J^ffg is the excess contribution to the hard sphere free energy over the ideal gas value. Finally, we remark 
that terms of higher order in F"^ representing dynamical corrections to the overdamped result, such as cj\^ (X, r), are 

functionals of (/iq"' (X, r), but cannot be expressed as functional derivatives of the free energy. They play a role only 
when the system is out of equilibrium, i.e. when currents are present. 

V. DISCUSSION 

In this paper we have derived the DDFT equation for the noise averaged density p(x, t) relative to a system of 
particles with inertia and driven by a stochastic heat-bath. The starting point is the many particle Kramers equation 
for the phase space distribution which corresponds to the Langevin equations of motion. Such a microscopic model 
describes fluids, where, due to the interactions with the solvent, the momentum and energy reach equilibrium on a 
faster scale than the slowly varying number density. In more detail, we employed the standard BBGKY hierarchy to 
derive an exact evolution equation for the one particle distribution in terms of the two-particle distribution function. 
In the case of a soft intermolecular potential we have applied the so-called Vlasov approximation to derive a closed 
equation for the distribution function. On the other hand, in the case of hard spheres, in order to truncate the BBGKY 
hierarchy, we employed the RET approximation, which is tantamount to factoring the two particle correlation into 
the product of two one-particle phase space distributions times the equal time equilibrium configurational two particle 
correlation function. The resulting equation is similar to the Enskog-Boltzmann equation for the one body phase-space 
distribution augmented by a Fokker-Planck term accounting for the heat-bath. 

In practical applications a realistic intermolecular potential can always be approximated by a hard-sphere potential 
and a soft potential. Hence, by a slight modification of the theory presented here, one can treat the hard-core repulsion 
with the method illustrated in section III.B, whereas the soft term is treated according to the method of section UFA, 
assuming that the processes induced by these forces proceed without interfering. 

The further approximation of this paper is to use an expansion in the inverse friction parameter to eliminate the 
modes associated with momentum and energy fluxes. In order to avoid the appearance of secular terms this procedure 
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is carried out using a multiple-time scale analysis. 

Interestingly, the method proposed yields, to leading order in the inverse friction parameter, an equation identical 
to that obtained in the DDF theory. Moreover, the merit of the present approach is the possibility to describe systems 
subject to externally imposed temperature gradients^^ or inelastic non conserving interactions^-' Finally, the theory 
allows to compute higher order corrections, whose next order to DDF theory represents the coupling between the 
density and the momentum current. Given the above, our result could be employed to study the dynamics of systems 
where inertial and structural effects result in a non trivial interplay, such as for systems undergoing Poiseuille flow or 
under shear. Candidate systems are colloidal suspensions and dense confined fluids. 

Finally, we wish to comment that the present theory can be useful for obtaining information about the relaxation 
of the structural properties, but does not allow to study hydrodynamic phenomena. In fact, Galilean invariance 
and conservation of momentum are explicitly violated by our model. The Langevin coupling dampens the absolute 
velocities due to the presence of a solvent at rest in the laboratory frame. However, the conservation of momentum 
could be enoforced by introducing a translationally invariant thermostat which damps the relative instead of the 
absolute momenta of particles. Actually, this strategy is at the basis of the Dissipative Particle Dynamics (DPD) — 
and could be embedded in the framework of the present approach. 

U.M.B.M. acknowledges the support of the Project COFIN-MIUR 2005, 2005027808. 



APPENDIX A: HERMITE POLYNOMIALS 

The j/-dimensional Hermite tensor polynomials, h'^^iy), introduced in equation (jl9p are tensorial generalizations 
of the standard Hermite polynomials^^. These are defined through the formula 

Explicitly the first polynomials read 

/i(°)(V) = 1 

The tensor polynomials are orthonormal with respect to the weight uj{V) defined by equation ([20|) . 

I dYu{V)h^l\\)h'^l\Y) = <5j' (Al) 

where is zero unless jjL — v and the subscript a is a permutation of the subscript /?. 
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